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Section 1 
ALGEBRA 


1.1. POWERS OR INDICES 
If n is a positive integer, 

x” == x XX XX... XxX acontinued product of n factors 
Otherwise, provided x + 0, 

@n=0,x° = 


(b) n = —p, where p is real and positive, x-? = ~ 


(c) n = p/q, where p, g are both integers, x?/4 = %/x? 
If x = 0, 0" = 0 for all values of n. 


Laws of Indices 


1.2, LOGARITHMS 


DEFINITION: If x = a°, then log, x = c. 


log (xy) = log x + log y 

log (x/y) = log x — log y 

log (x") = n log x 

log ("+/x) = (1/n) log x 

logiox = Pi X loge x = loge x/loge 10 = 0-4343 log. x 
(ec = 2718 Va 


1.3. IMPORTANT FACTORS 


x? — a? = (x + a) (x — a) 
x? + a® = (x + a) (x? F ax + a?) 
Key + 2? — SxyZ 
= (x + yt Zz)? + yy? + 2? — xy — yz — zx) 


1.4. QUADRATIC EQUATIONS 

If the quadratic equation ax? + bx + c = Ohas roots «, 8, then 
x? + (b/a) x + cla = (x — a) — B) = x? — (a + Bx + af 
Sum of roots = « + 6 = —bd/a 
Product of roots = af = c/a 
a, B = [—b + (6? — 4ac)]/2a 

Discriminant D = b? — 4ac 


D > 0, roots real and unequal 
D = 0, roots real and equal 
D < 0, roots complex and unequal (complex quantities, p. 23) 


1.5. REMAINDER THEOREM 

If nis a positive integer and ay, Gn-1,. « ., Gi, Ao are constants, 
Sf) = anx" + An-1 oe tax +a 

(x) is a polynomial of degree n. 
If f(x) is divided by (x — c), the remainder is f(c). © 
Thus (x — c) is a factor of f(x) if f(c) = 0. 

1.6. RATIO AND PROPORTION 


The ratio of a to b (written a:b) is denoted by the fraction a/b. 
a, b, c and d are in proportion if a/b = c/d 
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If a/b = c/d then 


a+b cxtd 
ad = bc - alee 
O50 cab d la+mb_ Ic+md 
a—b c—d pa+qb  pa+qd 


(where /, m, p, g are real numbers). 


1.7. VARIATION 


If k is a constant, 

Direct variation: if y oc x, then y = kx 

Inverse variation: if y oc 1/x, then y = k/x 

If z « x when y is constant, and if z oc y when x is constant 
then z o xy, ie. z = kxy when both x and y vary. 


1.8. SIMPLE SERIES 


An Arithmetic Progression is a series of numbers such that any 
term is derived from the previous term by adding to it a constant 
quantity, called the common difference (d). 

Let a be the first term of the series, then the nth term is 


l=a-+(n— 1)d 
and the sum of n terms is given by 


S, =a+(a+d)+(a+2d)+?..+(@+n-— 1d) 
= $n(2a +n — 1d) = tn(at+/ 
A Geometric Progression is a series of numbers such that any 
term is derived from the previous term by multiplying the latter 


by a constant quantity called the common ratio (r). The nth 
term is ar"~* where a is the first term. 


S,=atar+tar+...t+a'=al—r)di—nvn 
The sum to infinity = So =a/(l1—r) (|r| < 1) 


Harmonic Progression. A series of numbers a, b, c,. . . will be 
in Harmonic Progression if their inverses 1/a, 1/b, 1/c,. . . are 
in Arithmetic Progression so that 


l/c — 1/6 = 1/b — 1a 


Means of Two Numbers 
Arithmetic mean = A = }(a + 5) 
Geometric mean = G = 4/(ab) 
Harmonic mean = H = 2ab/(a + b) = G?/A 


Sums of Important Series 
Sum of the first n natural numbers = 4n(n + 1) 
Sum of the squares of the first n natural numbers 


= 4n(n + 1)(2n + 1) 
Sum of the cubes of the first m natural numbers = [4n(n + 1)]? 


1.9. COMBINATORIAL FORMULAE 

Any arrangement of n objects in a definite order is called a 

permutation. 

The number of permutations of n objects amongst themselves is 
nn —1)(n—2)...2.1=n! 

The number of permutations of n objects r at a time is 
"P, =n!/(n — r)! 

The number of combinations of n objects taken r at a time is the 


number of selections, each of r objects, without reference to 
their order within a selection, and this number is denoted by 


"Or (TF 
r 


(") = n! _, Hee)... (2 aD 
r)  rin—v! r! 


A 


The number of ways in which n objects may be divided into k 
groups containing respectively m, no, N3,. . .,M, Objects such 
thatn, #ngt...+t+m=nis 


n!/(ny!ne!. . .m,!) 


1.10. BINOMIAL AND MULTINOMIAL THEOREMS 


Binomial Theorem 


(x +a) = (7) + ("janie + (")a-208 tea. 


+: (")a-t¢ ee © (“Jarl 


(n is a positive integer.) (See also Binomial Series, p. 44.) 


The values of the binomial coefficients for different values of n 
are: 


Value of n Binomial coefficients 
1 1 1 
2 1 2 1 
3 1 3 3 1 
4 1 4 6 + 1 
5 1 5 10. 1008 1 


etc. etc. 


Note that the binomial coefficients for any value of may be 
obtained by summing consecutive pairs of coefficients in the 
preceding line. 


Multinomial Theorem 
(x1 + Xo t+... + Xx)” 


n! 
= >—--— XaM1xgN2. 6 XY 
Ny! na!l. . .m! 


kk 
where n = } n, and each n, takes integral values from 0 to a. 
i=1 


1.11. DETERMINANTS 


Second order determinant | a; 5, | has the value (a,b2 — azb1) 
a2 be 
Third order determinant | a, 6; c, | has the value 
Az be Co | aA, + O1B, +0,C, 
a3 bs C3 
where A, = ‘be C2 B, =] Cg Ae Ci = a2 be 
bs C3 C3 As a3 b3 


The solutions of the simultaneous equations a,x + biy = cy 


a2x + bey = C2 
are given by 
mca MMM Se 
Ci b, a Cy ay b, 
C2 be ag Ca a2 be 


and the solutions of a,x + biy + ciz = d, are given by 
a2x + bay + Coz = de 
a3x + bay + C3z = dz 


x ne y si > ae I 
d, b; C1 ay d, C1 a, by d, Qa, bx C12 
do be C2 a2 do C2 a2 be dz a2 be Co* 
ds bg cg a3 d3 C3 a3 bs dz a3 bz C3! 
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MENSURATION AND 
TRIGONOMETRY 


2.1. ANGLE 


One radian is the angle subtended at the centre of a circle by an 
arc equal in length to the radius r. 


27 radians = 360° 

7/180 radians = 1° 

1 radian = 180°/a7 = 57-3° 

Length of arc = s = ré (@ in radians) 


2.2. TRIGONOMETRICAL RATIOS 


Let 0 be the angle through which 
the radius OA would turn to 
reach the position OP, i.e. 6 = 
ZPON. 

(Convention: counter-clockwise 
rotation gives a positive angle 
and clockwise rotation gives a 
negative angle.) 


ZOPN = 90° 


ON and PN are directed lines in 
accordance with the convention 


for Cartesian coordinates so that 
ON is positive in the Ist and 4th quadrants 
and negative in the 2nd and 3rd 
NP is positive in the 1st and 2nd quadrants 
and negative in the 3rd and 4th 


DEFINITIONS 
sin 0 = NP/OP cos 6 = ON/OP tan 0 = NP/ON 


Algebraic signs of the above ratios in the four quadrants 
+ + “ 


sin cos tan 


Further Definitions and Identities 


cosec 6 = 1/sin 0 sec 0 = 1/cos 0 cot 0 = I/tan 0 
sin~ +x is the angle 0(— 42 < 6 < 32) such that sin 0 = x 
cos~*x is the angle 0(0 < 6 <z) such that cos 6 = x 
tan~*x is the angle 0(— $a < 0 < 4m) such that tan 6 = x 


Identities 
sin? 0 + cos? 0 = 1 1 + tan? 0 = sec? 0 
1 + cot? 0 = cosec? 0 tan 0 = sin 6/cos 0 


Complementary Angles A + B= 90° 
sin A = cos B cos A = sin B tan A = cot B 


Supplementary Angles A + B= 180° 
sn A=sinB cosA=-—cosB- tanA = —tanB 


2.3. MENSURATION 
Plane Areas 


Parallelogram Sides a, 6 and included angle «. Perpendicular 
distance between parallel sides of length a = h = b sin a 


Area = ah = ab sin « 


Trapezium Area = } (Sum of parallel sides) x Perpendicular 
distance between them 


Circle Radius r, diameter d = 2r. 
Circumference = 27r = awd Area = ar? = 4rd? 
If an arc subtends an angle 0 radians at the centre 


Area of sector = 4r?0 
Area of segment = 3r? (0 — sin @) 


Ellipse Semi-axes a, b; area = mab. 


Solids 


Notation: V = volume, C = curved surface area, S = total 
surface area. 


Solid Circular Cylinder 
h = height, R = radius 
V = wR*h C = 27Rh S = 27Rh + 27R? 


Hollow Circular Cylinder 
r = inner radius 
V = n(R? — r*)h 


Pyramid and Cone 
V = 3 (Area of base) x Perpendicular height 


Right Circular Cone 
h = height, R = radius of base, / = slant height = \/(R? + h?) 
= 47R*h C= Ri 
S = 7R? + 7R!/ 
Volume of frustum = 47h (R? + Rr + r?) 


where R, r are radii of the ends and his the height of the frustum. 
Sphere Radius R =$7R® C= 4nR? 

Volume of a zone = @7h [3(R? + r?) + h?] 
where R, r are radii of the ends and h is the height of the zone. 
2.4. SOLUTION OF TRIANGLES 


Triangle ABC: sides a, 6, c. Perpendicular height above the 
base a is h. Angle A is opposite side a, etc. R, r are the radii of 
the circumscribed and inscribed circles of triangle ABC. 


Semi-perimeter = 5 = }a+b+ 0c) 
Area = tah = fabsinC = \/[s(s — a)(s — b)(s — c)] = rs 
If C = 90°, c? = a? + b? (Pythagoras’ Theorem) 

Sine formula: a/sin A = b/sin B = c/sin C = 2R 


Cosine formulae: a? = b? + c? — 2bc cos A ete. 
a=bcosC+ccosB etc. 


sin 44 = aa Le = £= cos 4A = ‘| as © ©) ms @) 
nia = /2=2 in 0 


Area of a cyclic quadrilateral (sides a, b, c, d) 
V/[(s —a)(s—b)(s—c)(s—d)] wheres=3(a+b+c+d) 
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2.5. COMPOUND ANGLES 


sin (A + B) = sin A cos B + cos A sin B 
cos (A + B) = cos A cos B F sin A sin B 

tan A + tan B 
Gas ee Pe tan A ton B 
sin 2A = 2 sin A cos A 
cos 2A = cos?A — sin?A = 2 cos?dA — 1 =1 —2 sin?A 
sin 3A = 3 sin A — 4 sin® A 
cos 3A = 4cos?A — 3 cos A 
sin C + sin D = 2 sin 4(C + D) cos 4(C — D) 
sin C — sin D = 2 cos 4(C + D) sin #(C — D) 
cos C + cos D = 2 cos #(C + D) cos #(C — D) 
cos C — cos D = —2 sin 4(C + D) sin $(C — D) 


2.6. SPHERICAL TRIGONOMETRY 


Spherical Triangles 


The section of a sphere by a plane 
through its centre is called a great 
circle. A spherical triangle is the figure 
bounded by three great circles. 

The sides BC, CA, AB of the 
spherical triangle ABC are measured 
by the angles a, b, c, which they subtend 
at the centre O. The angles A, B, C are 
the angles between the planes of the 
great circles. 


Spherical excess = E=A+B+C—a 
Area of spherical triangle ABC = (A + B+ C — n)r? = Er? 
Sine formula: sin A/sin a = sin B/sin b = sin C/sin c 


Cosine formulae 
cosa =cosbcosc +sinbsinccosA ete. 


il 


Cotangent formulae 
cot b sinc = cos ccos A +sin Acot B etc. 


cos A + cos Bcos C 
cos @ = ——__—___—___—__ etc. 
sin B sin C 


cos $(a — b) 
cos 4(a + b) 
sin $(a — b) 

cos $(a + b) ying 
cos 4(A — B) 
cos $(A + B) 
sin 4(A — B) 
sin $(A + B) 
Ifa+b+c=2s 


shes Sik /Fe et (s — b) sin (s — c) 


sin b sin c 


sin s sin(s — a 
cos 4A = jr (s—4a) ) 
sin b sinc 


a (mea (s — b) sin (s — c) 


sin s sin (s — a) 


tan (A + B) = cot $C 
tan 3(A — B) = 
tan 3(a + 5b) = tan 4c 


tan #(a — b) = tan 3c 


IfA+B+C=2S 
— cos Scos (S — A) 


ew sin B sin C 
_ [cos (S — B) cos (S — C) 
siege J sin B sin C 


tan be a if — cos S cos (S — A) 
cos (S — B) cos (S — C) 
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Delambre’s Analogies 
sin $(A + B) cos $c = cos 4(a — b) cos $C 
sin 4(A — B) sin 4c = sin 4(a — b) cos $C 
cos $(A + B) cos $c = cos 4(a + b) sin $C 
cos (A — B) sin $c = sin 4(a + 5) sin $C 
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Section 3 


COORDINATE 
GEOMETRY 


3.1. DIVISION OF LINE IN GIVEN RATIO 

P (x, y) divides the straight line AB in ratio m:n, i.e. 
AP/PB = m/n 

With A (1, yi) and B (x2, ya), 


PO | + MXe _ ny + mye 
m+n m+n 


[For external division, the sign of m becomes negative.] 


3.2. EQUATIONS OF A STRAIGHT LINE 


Standard form m= tan 0 = gradient or slope, c = intercept 
on y axis 
y=mx+e 


Perpendicular form p= length of perpendicular from the 
origin O to the line and « is the inclination of the perpendicular 
to Ox 


x COS + ysin « =p 
Intercept form a, b are the intercepts on Ox, Oy respectively 


42 
Py ' 
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General form 
Ax + By +C=0 
y—yi =m(x—x,) [Line of slope m passing through 
(x1, ¥1)] 


(x — x1) [Line passing through (x1, yi) 


yak ) 
and (x2, ya)] 


*g™" X1 


yn =( 


Angle between two lines y = mx + c, and y = mox + cz is 


-, ”™— Me 
tan ~? ——___>_ 
] + My, Moz 
Two lines are perpendicular if mymz = —1. 


Laws from Plotted Results. Form of law modified to give 
straight line graph. 


Law Other Quantity Quantity Slope Inter- 
suspected form of to be to be gives: cept 
law plotted _ plotted gives: 
along along 
horizon- _ vertical 
tal 
y=ax? +b ——__—? y a b 
xy = ax 
+ by yobylxta yix y b a 
iy = ax" logioy = logi0x logioy n logioa 
nlog x 
+logioa 
y = ae™ logioy = x logioy 0:4343k logioa 
0-4343kx 
+ logioa 
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3.3. CIRCLE 
(i) x? + y? =a? Circle centre origin, radius a 
(ii) (x — h)? +(y — k)? =a? $Circle centre (h, k), radius a 


(iii) x? + y? + 2¢x + 2fy +c =0 General equation of 
circle 


Tangent at (x, y1) on circle (i): 

XX1 + yy = a’ 

Tangent at (x1, y:) on circle (ii): 

(x — A(x. — A) + (y —kK)(n — bk) = 2? 
Tangent at (x1, yi) on circle (iii): 

xXXy + yyy te(xt+xm)+fly+y) +c =0 


3.4. CONIC SECTIONS 


A conic section is the locus of a point which moves in a plane 
such that the ratio e (called the eccentricity) of its distance from 
a fixed point (the focus) to its distance from a fixed straight line 
(the directrix) is constant. 

The locus is an ellipse, a parabola or a hyperbola according 
as 


e<=>1 


Polar equation of a conie referred to its focus where / = semi- 
latus rectum 


Dae aoe 
; 


(a) Ellipse 


Ellipse, centre origin, semi-axes, a, b, e < 1 


b? = qa? (1 — e?) [= 5?/a 
16 


Tangent to above ellipse at (x1, y1) on it 


XX V4 
he 


(b) Parabola 
y* =4ax Parabola, origin at vertex, Ox axis of symmetry 
yyi = 2a(x + x1) Tangent at (x1, y1) on the parabola 


y=mx + < Tangent of slope m for any value of m 
I 


Parabolic arc between two points at the same level. Span J, 
height s at the mid-point referred to origin at one end: 


4 
y= axl — x) 
Area of segment of parabola = § s/ 
(c) Hyperbola 


5-7 = 1. Hyperbola, centre origin, semi-axes a, b 
> — = 1 Tangent at (x1, y:) on the hyperbola 


~ + $ = Equations of the asymptotes 


b? =a? (e?7—1) (e>1) 


xy =a’ Rectangular hyperbola, asymptotes are the coordinate 
axes 


YiX + X1y = 2a? Tangent at (x1, yi) on rectangular hyperbola 


(d) Focal Distances 


The ellipse and the hyperbola have two foci and two directrices 
whilst the parabola has only one focus and one directrix. A 
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focal distance is the distance of a point on the conic from a 
focus. 


Ellipse: the two focal distances are a + ex. Sum = 2a. 
Hyperbola: the two focal distances are ex + a. Difference = 2a. 


3.5. GEOMETRY OF SOLIDS 


(a) Point 


Let (x, y, z) and (r, 0, 6) be the Cartesian and Polar coordinates 
of a point P, 


x=rsin@cos¢ y=rsin@sin ¢ z= rcos 0 


where r = OP = +/(x? + y? + 2?) 
6 =co-latitude and cos 0 = z/./(x? + y? + 2?) 
@ = longitude and tan ¢ = y/x 


(6) Straight Line 
Distance between two points P (x1, yi, 21) and Q (x2, Yas Ze) 
d = V[(x%1 — x2)? + (v1 — yo)? + (21 — 22)?] 


The Direction Cosines of a straight line PQ are the cosines of 
the angles between PQ and Ox, Oy, Oz and are denoted res- 
pectively by /, m, n. 


fe NEM Yom Vi Zg— 21 


n= r= 


d d d 
f2+nm?+n2=1 


The angle 6 between two straight lines whose direction cosines 
are (/;, m1, m1) and (lz, me, na) is given by | 


cos 9 = Lylg.+ MyM + NyN2 
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Equation of the straight line passing through the point (a, b, c) 
having direction cosines (/, m, n) is 


(x, y, Z) = 0 represents the equation of a surface 
(c) Plane—linear in x, y, z 
General equation of a plane 
ax + by +cz+d=0 
Intercept form: plane cuts the axes at distances a, b, c from O 
Cone oo | 
aes « 
Perpendicular form: perpendicular from O to the plane has a 
length / and direction cosines (/, m, 1) 
Ix + my + nz =p 
Length of perpendicular from (x1, yi, 21) to the above plane 


Ix, + my1 + NZ, 7 
(d) Surfaces of the Second Degree 


Surfaces of revolution 


x? + y® =a? Right circular cylinder, radius a 
x? + y® = a®z° Right circular cone, vertex origin, axis Oz 
x® + y? 4+ 2% =a? Sphere, centre O, radius a 


ry ae ae did ls : 
* + a + , = 1 Ellipsoid, semi-axes a, b, c 
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Section 4 
NON-DECIMAL ARITHMETIC 


4.1. ARITHMETIC BASES 


Decimal Notation: numbers are expressed as multiples of powers 
of 10. 


10 is called the base. The symbols used are 1, 2,. . ., 9. 
293 means 2 x 107 +9 x 10+ 3 


Any other positive integer may be used as the base of a number 
system, e.g. in octal arithmetic, 537 means 5 x 82 + 3 x 8 + 7. 
Decimal equivalent: 351 

More generally, in arithmetic base r, the number represented by 


Qn An-1- « «4281, Qo 
where each of the digits ao, @i, da,. . +, @n-1) Qn < r(Gn # 0), is 
Ant” + Anaar™™* +... Hb aar? + ar + ao 


Binary arithmetic (base 2), using the two digits 0, 1 only, has 
become very important because it is widely used in digital 
computers. 


4.2. CHANGE OF BASE 


N is an integral number expressed in any base. It is required to 
express N in base r. 
Let N = ayr” + dn-r™ 1 +... bart ao. 
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Divide N by r so that N/r = agr*~* + Qp-ar™-? +... + a1: 
remainder do. 

Divide the quotient by r. The next remainder is a, and so on. 

The successive remainders are the required digits ao, a1,. . ., Qn. 

For example, by successive divisions by 2, decimal 247 becomes 

11110111 in binary. Binary numbers are most easily converted 

to decimal by reversing the process 


Looe Oo OR OT wt 
Sigs. 30 61 125) 247 


Multiply the left-hand digit by 2, add the next binary digit, and 
continue this process. 


4.3. FRACTIONS 
2 6 7 


In decimal arithmetic, 9-367 9+ — + — + —. 
n decimal arithmetic, 9 means 9 + 10 * 102 * 103 
Similarly in binary arithmetic, 1.101 means 1 + 5 + — + 5 
In arithmetic base r, the fraction 0. 6; b2b3. . . where 5, bo, ba, 

. . are integers, each less than r, of which one or more may be 
zero, is equivalent to 


Oi 5 Be. Be 
, + p2 + Pc) + eee 
To express any fraction F in the scale of r, b;, bz, bs, . . . must 
be found such that 
pat y bey bey 
‘RAT a r 
Multiplying byr Fr=6, + 22 +8 Pie 0 


so that b, is the integral part of Fr, the fractional part F, being 
given by 


and bz is the integral part of Fir and so on. 
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For example, to express 11/16 (decimal) as a binary fraction: 


X2=1+h #xX2=042 EX2=144 
Pee Lie 


11/16 (decimal) = 0.1011 (binary) 
4.4. THEOREM 


In the scale of r, the sum of the digits of any integer when 
divided by (r — 1) gives the same remainder as when the 
integer itself is divided by (r — 1). 


Corollary: in the decimal scale (r = 10), a number is divisible 
by 9 if the sum of its digits is divisible by 9. 
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Section 5 
COMPLEX NUMBERS 


5.1. COMPLEX NUMBERS 


If a and b are real numbers and i = «/—1 (i.e. i? = —1), ib is 
known as an imaginary number and z =a + ib is called a 
complex number and Z = a — ib is called the conjugate of z. 

If z, =a+ib and z,=c+ id, their sum, difference and 


product are 
%m+2=@tc)t+iO+td). 
212g = (ac — bd) + i(ad + bc) 


z2Z = q? + 5? 

a _atib_(at ibe—id) 

Z2 ¢ctid (c+id)(c — id) P(oc, y) 

_ (ac + bd) + i (be — ad) : Mi 
c? + d? 

If z, = Zo, thena=c,b = d. | 8 a 


O 


5.2. ARGAND DIAGRAM Argand diagram 


There are two geometrical ie "napagegamm of the ea 
number z = x + iy: 

(a) the point P(x, y) (6) the vector OP 
| 2=x+iy=r(cos 0 + isin 4) 
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DEFINITIONS 
(a) Absolute Value or Modulus of z = |[z| = |x + iy| 
= length of OP = r = (x? + y’)1/2 
(6) Argument or Amplitude of z = arg z or amp z = 9 
= tan? (y/x) where —7 <0 < a. 

i OR Ra oe.) ee 
of oe La ar » Ree 


a 6? 6 
= (1 ae ae ° ) +i(0-T4+o-. ° .) 
= cos 0 + isin 0 
Therefore e~‘® =cos@— isin @ so that 
cos 0 = 3(e? + e- 9) = cosh i 
isin 6 = $(e'? — e- %) = sinh i0 
Therefore i tan 6 = tanh i0 
Also z=re® and Z=re~‘ 
21Z2 = ri(cos 6, + isin 6) . re(cos 02 + i sin 6.) 
= 7 e'91 , roe! 
= ryrze'1 +92) = ryra(cos 0; + 02 + isin 0, + 82) 
arg 2172 = 0, + 02 = arg 7, + arg Ze 
[Z1Z2| = rire = |2,| . [Ze] 


5.3. DE MOIVRE’S THEOREM 


{r(cos 0 + isin 6)}" = r"(cos nO + isin nO) for all values of n 


nth roots of z 
z=x+iy =r(cos 6 + isin 0) 
= r{cos (0 + 2km) + isin (0 + 2km)} = ret(O+2km) 
where kK = 0, 1,2,.. ., (2 — 1). 
Rf/z=/(x +iy) = rin (eos | ae + isin t ra a) 


n 
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nth roots of unity 
1=cosO0+isinO _ therefore 


"4/1 = cos (==) + isin (=) 


e. 
Cube roots of unity (7 = 3) 
*yi mec (2) +n (2) 


The three roots are 1, —$ + i/3/2. 


[A =0, 1, 2, 


(kK = 0, 1, 2) 


..4(4— 1)] 
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Section 6 
MODERN ALGEBRA 


6.1. SET 


A set is a well-defined collection of objects. ‘“Well-defined” 
means that it may readily be decided whether an object is a 
member of the set or not, e.g. 


(i) The set of all counties in the U.K. The elements of this set 
may be listed: {Aberdeenshire, Anglesey, . .. Yorkshire}. 
(ii) The set of all positive even integers less than 16: {2, 4, 6, 
«a ee 


A set may contain an infinite number of elements. 


6.2. NOTATION 
x € A means x is an element of the set A. 
(x ¢ A is the negative statement.) 


A = {x1, X2,. . .» Xn} Means x1, X2,. . ., Xn are all the elements 
of set A. 


A = {x, P} means A is set of elements x satisfying a property P. 
43 means “‘there exist(s)”’ 
For two statements p and q, 

Pp > q means “p implies g” or “if p then q”’. 


p <> q means p and g are equivalent statements or “‘p if and 
only if q’’. 
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A = B means two sets A and B are equal, i.e. A and B contain 
precisely the same elements or x € A > xe B. 


A ~ B means sets A and B are equivalent, i.e. there is one- 
to-one correspondence between their elements. 


A © Bmeans A is a sub-set of B, ie. every element of A is 
also an element of B. 


Ac B<+(xEA> xEB) 
A=B<+Ac BandBCA 


6.3. OPERATIONS ON SETS 


A U B denotes the union of two sets A and B. It is the set of 
those elements which belong to either A or B or both 
A and B. 
1ne.xE€AorxEeEBsoxEAVUB 

Ac B denotes the intersection of two sets A and B. It is the set 
of those elements which are members of both 4 and B. 
ie.xE€AandxEBsxEANB 


The universal set V is the set of all possible objects in a given 
problem. The empty or null set ¢ is the set having no elements. 
If two sets A and B have no elements in common, 4 Nn B = ¢., 


6.4. VENN DIAGRAMS 


Let the universal set V be represented S' 
by a rectangle and let any subset S$ V 
of V be represented by.a circle 
(shaded) drawn inside the rectangle. 
The set S’ (unshaded) of all elements 
not in Sis called the complement of 
Sand Su S’ = V. 


Aj Bis represented by the doubly 
cross-hatched area. 


A U B is represented by the whole 
shaded area 
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The following important results can easily be verified by Venn 

diagrams: 

QM (4A OB) ANC=AnN(BAC BAW : 
(AUB)UC=AU(BU a Associative Laws 

Gi) AU(BNQ=(AVYVBN(AVUOC) 
AN(BUC)=(ANB)UANC) 

(ill) a ye ee if . 2 m De Morgan’s Laws 


Distributive Laws 
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Section 7 
MATRICES 


7.1. THE ORDER OF A MATRIX 


A matrix is an ordered array of numbers or other elements. 
A matrix having m rows and n columns is said to be of order 
mxX Nn. 


Qi1 G12 Aig 2\ is a column matrix; 
Q21 A222 ao23 1} order 3 be 
Q31 A322 Ass 3 


G41 442 Gag 
is a4 xX 3 matrix. 
(—2 1 O 7)is arow matrix; order 1 x 4. 


A matrix having 1 rows and n columns is called a square matrix 
of order m. Matrices will be denoted by heavy capital letters 
ye, By 6s > 


7.2. THE MULTIPLICATION OF MATRICES 


Suppose A is a 3 x 5 matrix and Ba 5 x 4 matrix. Then, for 
example, the element c23 of the product A x B is obtained by 
the multiplication of the elements of the second row of A by the 
corresponding elements of the third column of B, adding the 
products so that 


Cos = Aaibig + Agabe, + Aagbg33 + da4bag + do5b53 


5 
= 2) A21Dx3 
k=1 
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In general 


5 
Cy D> AnD 
k=1 


Compatibility 
The above product rule is applicable only if the number of 
columns of A equals the number of rows of B. Matrices which 
satisfy this criterion are said to be compatible for multiplication ; 
otherwise, their product is not defined. 

In general, a matrix A (m xX n) is compatible with a matrix 
B (a X p) and A X Bisa matrix C (m xX p). 

A matrix A (m xX n) is compatible with a matrix B (x x m) 
for multiplication in either order but A. B 4 B. A in general. 


7.3. UNIT MATRIX: NULL OR ZERO MATRIX 
( tO 4 =( A = (2 sib 4) 
0 1/\c d ce @ ec. dye 1 
If any (2 X 2) matrix is pre- or post-multiplied by I = . Os 


it is unchanged. 
I is called the unit matrix of order 2. 


1 0 0 
(0 1 0) is a unit matrix of order 3, and so on to higher 
me 814 


orders. 
A zero matrix has all its elements zero. 


7.4. EQUALITY; ADDITION; MULTIPLICATION BY 
A SCALAR 


(a) Two matrices are equal if and only if they are of the same 
order and their corresponding elements are equal. 
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(6) The sum of two matrices of the same order is defined to be 
a matrix whose elements are the sums of the corresponding 
elements of these two matrices. 

(c) The product of a matrix by a number (scalar) k is defined to 
be the matrix derived from the first by multiplying each of its 
elements by k, e.g. 


Lee O _ {ka kb 
ra =(¢ 7) then kA = (12 ta) 


7.5. DETERMINANT OF A SQUARE MATRIX 


Matrix A (n X n) Determinant of matrix A = |A| 
11 Aig. ..-Ain 411 Gyo... .Ain 
a21 Q22 aan G21 A22 Aan 
Ani Ane Qnn Qn1 an2 Ann 


A rule assigns a value to the determinant |A| of order n. 


Let a,,; be the square matrix of order (n — 1) obtained by 
deleting the ith row and the jth column of A. The determinant 
|a,s| is called the minor of ay in |A|. The cofactor Ay of ay is 
defined by 


Ay = (—1)'*! X value of [a,,| 
The value assigned to the determinant |A| is 
n 
JA] = 2% ay Ay 
j= 
(This formula will give the same value of |A| irrespective of the 
particular row in terms of which it is evaluated.) 


If the products obtained by multiplying the elements of any 
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row by the corresponding cofactors of another row are added, 
the result is zero. 


nr 
ay, Ay =O (i #f) 
k=1 


ae fe son 


(2x 2) G21 422 
a1 a 
|A] = sn fic = 411411 + 12412 = 11422 + Ai2 (—aa1) 
21 22 


Q4i1 Qi2 Aig 
A=1[do1 Qe2 ag 


(3*%3) \@g1 Ase 33 
Qi1 Qi2 Aig 
JA] =| @e1 de2 Gaz 


Q31 QAsq 433 
= 411Aq1 + 12A12 + 413413 


Q22 Aas 
QA32 a33 


G21 422 
Q31 32 


Q21 Qo2z3 
Q431 433 


= a1 


A matrix whose determinant is zero is a singular matrix. 


One of the rules for multiplying two determinants of the same 
order is precisely the rule (Sec. 7.2) for multiplying two square 
matrices of the same order. Hence if A, B are square matrices 
of the same order 


|AB| = |A]. |B] and |AB| = |BA| 
even though, in general, AB ¢ BA. 
7.6. INVERSE MATRIX 


The inverse of a square matrix A is A~+ of the same order such 
that 


A-1A = AA“? =I 
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eo ie b t= 7 d 4, 
Ira =(¢ 4)» then A =a he if 


A singular matrix has no inverse. 


7.7. MATRIX SOLUTION OF LINEAR EQUATIONS 
Equations ax+by=p and cx+dy=q 


Matrix form: (¢ 4 (*) aN (?) 
ated aie a q 
Solution: (* ghd od 4 Kt v) 
Jy [A] \9q — <p 
7.8. THE TRANSPOSE OF A MATRIX 


If the rows and columns of matrix A (m x n) are interchanged, 
the resulting matrix A’ (7 x m) is said to be the transpose of A. 


A square matrix A is symmetric if A = A’ and skew symmetric if 
A’ = —A. 

If A = A’ (symmetric), the elements a,; of A are such that 

ayy = Ay. 

If A’ = —A (skew-symmetric) a,; = —ay, (i 4), whilst a,, = 0. 


7.9. ADJUGATE OR ADJOINT MATRIX 


If A is a square matrix, the adjugate matrix adj A is formed by 
replacing each element of A by its cofactor in |A| and transposing, 
so that 


adj A = (Ajj)’ = (Ay) 
where (A,,) denotes a matrix of order n having elements Ay,. 


Q11 Gig... Ain Aix Aor... Ani 
Q21 A2e2 dan Ayo ~Aoe Ane 
then 
fA= F* : 
’ . adj A = 
QAn1 Ane Ann Ain Aon Ann 
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7.10. INVERSE OF A SQUARE MATRIX 


If A is a non-singular matrix of order n, 


eee 
A-+? =— adj A 
fal 
7.11. SOLUTION OF LINEAR EQUATIONS 
(@) QyX%1 + QyaXo t+... Fanti =O, G=1,2,...,”) 
are n simultaneous linear equations in m unknowns x. 
Matrix form AX = B 


where A is a square matrix of the coefficients a,; and X, B are 
column matrices of the x; and 5. 


If A is non-singular, XK = A~*B where A~? is given in 7.10. 


(b) Method of Triangular Decomposition illustrated for three 
unknowns. 


Matrix form AX =B (Ais the 3 xX 3 matrix of aj) 


1 0 0 UW41 U2 U3 
Find V = ( 1 0} and U = (0 Uz2 i 
0 O 


Uzi Us 1 Ugg 
such that VU = A. 
From the equations, Uji = Q11, Uig = Ala, . - +» Vai U1 = Gai, 
J *9 etc. 


Let VY = B(Y is a3 X 1 matrix) whence y; = 5; and 
VaiVi + yo = ba V3i V1 + UVs2Y2 + Ys = bs 


from which Y is found. 

AX=B VUX=VY UX =Y 
whence W11%1 + Ui2X2 + UisX3 = yi and 
Uz2X2 + U23X3 = V2 UssX3 = Ya 

from which X is found. 
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Section 8 
VECTORS 


8.1. COLUMN VECTOR 
— 
Position vector of P(a, b) = OP = ry. 


r, may be represented by a column matrix or column vector fet ; 
Addition of Two Vectors. Let rz = (5) , then 


mtr = (5) “++ (4) r= ( / id 4) (Parallelogram Law) 


y R(a+c, b+d) 


Base Vectors 


i = ( a and j = i represent unit vectors along Ox, Oy 


respectively. 
r= + =a(4) +6(9) = ai + bj 
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8.2. LINEAR TRANSFORMATIONS IN A PLANE 


A column vector (3) is transformed into another column vector 


(7) on pre-multiplying the former by a 2 x 2 matrix A 


(transformation matrix). 
Pere ai 
y y c d} \y 


ie. r’ = Ar and x’ = ax + by, y’ = cx + dy 


The base vectors i, j transform into i’, j’ such that 


wi=(2 4) (b)=(2)=# amt w= (0) 7 


i.e. i’, j/ are the columns of A. Thus if the images of the points 
(1, 0) and (0, 1), under a transformation, are known, the 
transformation is known. 


8.3. AFFINE TRANSFORMATIONS 


x’\ _ fax+by\ _ a er °, 
Siting) 31) +7 (a) aay 
The point (x, 0) transforms into (ax, cx), i.e. x (a, c). 
The point (0, y) transforms into (by, dy), i.e. y(b, d) 


A grid of unit squares in the (x, y) plane will transform into a 
grid of parallelograms. The transformation is such that parallel 
lines remain parallel and the ratios of lengths in the same 
direction are unchanged. It is called an affine transformation. 

Exceptionally, if the matrix A is singular, i.e. if |A| = 0, then 
alc = b/dor a = kb, c = kdso that x’ = (b/d)y’ and all points 
in the (x, y) plane are transformed into points on one straight 
line. 
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8.4. BASIC TRANSFORMATIONS 

Type Matrix Algebraic 
Equivalent 

Reflection in Ox 

Reflection on Oy 

Reflection in line y = x 


Reflection in line y = —x 


( 
( 
( 
( 
Quarter turn about origin i ey 
( 
( 
( 


Bo x 

Half turn about origin wei 4 ia Mase 

" * as a | x= y 
Three-quarters turn about origin hy ‘) yo 
Enlargement . y 7 5 

, 1 k x= x aad ky 
Shearing parallel to Ox ( 0 A ps) 
Rotation through angle « (ee % sin 

sin «% COS % 


x’ =x cosa — ysina 
y =xsina+ycosa 


Reflection in line y = x tan « ( Cos2% sin 2% 
sin 2a —cos 2a 
4 


x" = x cos 2% + y sin 2a 
y’ = x sin 2x — y cos 2a 
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8.5. EIGENVECTORS 


Consider the transformation 


OM) iG) 
y A ons 88 
where the transformation matrix A is non-singular. Vectors 


(7) whose directions are unchanged by the transformation are 


called eigenvectors, i.e. ( 0 - eo 
Ss" y ky 
We require solutions of ig b (3 ) oi k(* 
i.e. of (A — «D)(*) = 0 
The eigenvalues of A are the values of k obtained from the 
determinantal equation 


|A — Al|=0 
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Section 9 
DIFFERENTIAL CALCULUS 


9.1, LIMITS 


A function f(x) is said to approach the limit L as x approaches c, 

ic. lim f(x) = L if, given any positive number «, there is a 
Yaa #4 

positive number 6 such that 


If@ —L|<e when 0 <|x —c| <6 
_ Some important limits 


sin 0 __ tan 0 
lim ——. =_ | = | 
0-0 0 0-0 0 
lim sciatic l 
r->0 x 
i ee x a3 
wma x= 7 | 
lim (1 ++) = € lim (1 +2] == e* 
n— 00 n n— 0 n 


9.2. DIFFERENTIATION 
Let y = f(x) and let dy, dx be corresponding inorements in y,x 
then the differential coefficient or derivative of y with a tox 
is given by a Hae 3 
4 tim 8 a tim £242) — fe) 
dx Seo Ok . aero Ox 
_ = gradient of the curve y = f(x) at the point (x, y) on it. 
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Other symbols for this differential coefficient are y’, Dy, and 
f°@). 


Rules for Differentiation 


Let u, v, w, etc. be functions of x. 


y = cu (c = const) J ach 
OR dy du , dv , dw 
D iE he at a de ae ee ge oh 2 
me dy dv du 
y = w ao ps ag ae 
wie Y = 4 (oS — v2) 
+78 a a dx 
a dy. es 1 dv ee 
sepia Ye ig a op aa w dx 
— 7" dy _ n-1 du 
Se eee 
ay ay du df 4d 
y=fw Ee A a Mie 


where u = ¢$(x) 


If y and x are both functions of , 2 ws 2 - 
It x,y are distances and ¢ is time, sad 2 (often denoted by 


x,y) are velocities. 
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Table of Important Differential Coefficients 


y dy|dx 
x” nx” -1 
eks kekz 
sin (ax + b) a cos (ax + b) 
cos (ax + b) — asin (ax + 5) 
log. (ax + b) al(ax + b) 
tan x sec? x 
cot x — cosec? x 
sin~1 x/a 1/./(a? — x?) 
cos~? x/a — 1/\/(a? — x?) 
tan~? x/a al(x? + a?) 
sinh x = #(e*’ —e-*) coshx = ke" + e~*) 
cosh x sinh x 
tanh x aha sech? x = 1/cosh? x 


e* ob e~ 
9.3. NEWTON’S METHOD OF APPROXIMATION 


If xo is an approximation to a root of f(x) =0, a better 
approximation x, is given by 


so that X1 = Xo — f(xo0)/f (0). 
9.4. SUCCESSIVE DIFFERENTIATION 


If dy/dx is differentiated with respect to x, the second deriva- 
tive d*y/dx? is obtained. Successive differentiations produce 
derivatives of higher orders. 
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Notation 

=o. 2 (2) 
Ka Aes Se dx \dx] — dx* dx 
d (o2) 0) 
ax: \de 1 a 
nth derivative = £2 = Dry = fx) 


[f'@)] =f’@) 


Leibniz’ Formula for the nth derivative of a product. 


Dw) = D"u.v + (7) De-tu, Dov + (5) De-tw, D2vu+... 
n de n ie 
+ (7) >" "EPO ae fea be. D* ty +u. D*v 


9.5. MAXIMUM AND MINIMUM VALUES 
If dy/dx = 0 at x = a, then 

y = f() has a minimum at x = a if d*y/dx? > O atx =a 

y = f(®) has a maximum at x = a if d?y/dx? << Oatx =a 
If d?y/dx? = 0 at x = a, the graph of y = f(x) hasa point of 
inflection at x = a. 


9.6. CURVATURE 
y = f(x) graph APB 


PT, PN are tangent and normal 
at P respectively. 
dy/dx = tan yp 

= slope of tangent 


EES 
Subtangent QR aya 


= yly’ 


Subnormal QN = y 2 = yy 
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Length of tangent = PR = y(1 + y)*?/y’ 

Length of normal = PN = y(1 + y”)?/? 

Let s be the length of the arc AP measured from a fixed point A. 
Curvature = « = a =y" (1 + y?) 9? 
Radius of curvature = 1/x 


If the curve is nearly horizontal at all points as in the case of 
horizontal beams, dy/dx is small compared with unity and the 
curvature is d*y/dx? approximately. 


9.7, PARTIAL DIFFERENTIATION 


. se OF cal ch OE tah Beco 
Ma = [0 Y) 0s ee a co 
dx, dy, dz are called differentials. 


ff hE ale 
dt oxdt oydt 
If u = f (x1, X2, X3,~. » +, Xn), then 


If x, y, z are all functions of ¢, 


tig — Et gee al SO ay, aah Tt. ee dX 
Ox, OXe2 OX3 OXn 
Implicit Functions If y is an implicit function of x as given by 
f(x, y) = 0, 
of of dy of /of 
= J oe = ep a ee 
oe a ee 
9.8. INFINITE SERIES 


A series may contain either a finite or an infinite number of 
terms. Let S, denote the sum of n terms of a series. The series 
is said to be convergent or divergent according as the lim Sy, is 
a finite quantity or is + ©. idl’ 
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Binomial Series With reference to 1.10, if n is a positive integer, 


ee n n n 2 n .T 
ia a =(5) + ({)2* (3) et Ne Fs 
a, (ae 
<1 txt n(n — Dx? +. . t+ 


< n(n — 1). ~. (n—rtdDxt+...4+2x 


which is’a finite series. 
If n is not a positive integer 


(+x = 1 nx + Sala — Dx? +... + 


Sn(n — 1). --(2—r+i1)x' +. . .ad inf. 


provided |x| <1, ie. the infinite series converges to (1 + x)" 
provided |x| < 1, otherwise it is divergent. 


Approximations. If x, y, z, u, v are small compared with 1, 
(1 +x)"+1+nx = foralln 
e.g. /(1 + x) = 1+ 4x 


d+ od + yl + z) 


| oi Se 
(+00 +0) rues ee oe 


Exponential Series 


ef = 1 + (kx) +5 (kay 4 = (kay? . 7 (ea! aa ad inf. 
SP pep Se) ee 
e= 7 31 a ay = arity. ane 
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Trigonometrical Series 


_ 1 5 ; 
sin x xp toes, . .ad inf. 
1 . 
8 x = 1 ee +7! +...ad inf. 
“nx ech eee i oe +...ad inf. 
15 315 


Maclaurin’s Theorem 


f) ={O +O + 527'O+. +527" OH... 


Taylor’s Theorem 
Single Independent Variable 


fe +h) = f(x) + hf’) + =f) maka! + <hr f(x) 
ei 


Two Independent Variables 


LAT YT Re (az +4) 
2 OF o*f Pea 
+3 (we 5 + hk ae +k aa) + 


9.9. MAXIMUM AND MINIMUM VALUES 
By Taylor’s Theorem, neccesary conditions for a maximum or 
a minimum of f(x, y) are 0f/0x = of/oy = 0. 


2 ofr 2 
Provided that oa: ss > (=) , then f(x, y) is a maximum 


or minimum according as 6°//0x? and 6?f/@y* are both negative 
or both positive. 
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Section 10 
INTEGRAL CALCULUS 


10.1. STANDARD INTEGRALS 


f@ Sf dx 

Di x"*1/(n + 1) except when n = —1 
1/x logex 

i/(ax + 6) = log. (ax + ) 

p 1 

oe a®*/k log. a 

sin (ax + b) o * cos (ax + b) 


cos (ax + b) = sin (ax + b) 


tan (ax + b) log sec (ax + b) 


sec? x tan x 
cosec? x — cot x 
sec x tan x Sec Xx 
cosec x cot x — cosec x 
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f@) Sf (x)dx 


1/./(a? — x”) — sin=+ (x/a) 
1/VG2? +42) log : [x + V(x? + a)] = sinh-*(x/a) 


1/V@2 —a%) log ; [x + V(x? — a)] = cosh-*(x/a) 
1/(x? + a?) : tan~+(x/a) 


1/(x? — a?) - log (==4) (x > a) 


2a x + 
eee 1 atx 
I/(a x?) 5, [08 (<2 (x <a) 
sinh x cosh x 
cosh x sinh x 
sech? x tanh x 
flO) log f(x) 
fOOl/fO) 2 f(x) 
Integration by Parts 
cle OC 
fuZ ax =w foax 


Hence f log (x) dx = x log x — x 
10.2. REDUCTION FORMULAE 


a sf cos"~2 @ do 
n 


| cos 6 d0 = = sin 6 cos"=1 8 + 


n—1 


{ sin 6 a0 = — cos 0 sin*=* 0 + [ sinr-2 0 ao 
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f sin™ 6 cos" 6 dé 


1 — 1 
= — sin™~+ @ cos"*1 @ 
m+n ae 


at . sin"** 9 cos*-* 6 + = 
m+n 


— | sin”™~? 0 cos" 6 dé 
i 


=| sin” 0 cos"~? 6 dé 


| tan" 0 d@ = 


; ; ‘ane 7] mr ane 


| cot 6 ao =— . cot? — [ cot? 6 do 


10.3. iis INTEGRALS 
a /2 


27n/p ‘ 
\ noes sin pt dt = | a cos ptdt = 0 ninteger 


a /2p 
‘ ” sin pt dt = 1 cos pt dt = 1/p 


i sin? pt dt = shine cos? pt dt = 7/4p 
lig cos” 0 d0 =" Heh ine cos"~? 6 d@ 


Ee wo 2 | sin"~? 0 do 


/2 ine 
1 sin” 0 cos" 0.d9 = ™—! * sin™-2 9 cos" 0 do 
0 m-+-n vo 
nian /2 ; 
eon” sin™ 0 cos"~*,0.d0 
m+n Jo 


If p and q are integers, p # q, 


i : yO 
> Ssinpe sing? dt = 9 Sin pt cos qt dt 
Qa 
=}, cospt cos qt dt = 
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If p = q these become respectively. 
ng sin? pt dt = 7 lio § sin 2pt dr = 0 


on 
| cos? pt dt = 7 
0 


10.4. AREA BOUNDED BY A CURVE; LENGTH OF 


CURVE y 
ds y = f(x) 
Cartesian Coordinates S 
Area A = i foods 
1 
aS ie | “2 / + y'Ddx 
«nga %1 J O xy doc Lo ad 
Polar Coordinates nn 82) 
Area A = | 2 442 dO , 
se 01 i (4,9 
0 
Arcs = | i /[r? + (dr/d6)?]d0 
Z. 
= [ vil + olayrtar 
1 
0) 


10.5. VOLUMES AND SURFACES OF REVOLUTION 


When revolved about Ox, the curve y =f (x) generates a 
volume and a surface of revolution. 


ra 
Volume = | my*dx 
zy 


Area = | 2r yds = \ c 2ryV(1 + y)dx 


10.6. SIMPSON’S RULE FOR AREA AND VOLUME 


Area Under a Curve Let the area be divided up into an even 

number n of equally wide strips (width h) by ordinates y,, 

V29- + +> Vnt1 then 

A= th [(y1 + Yn+1) + 2(ys Tyst Yet... ) 
+4yotytyet...] 

Volume of a Solid Let the solid be divided up into an even 

number n of slices by (n + 1) parallel planes at equal distances 

h apart and let Ai, Ao,. . . Any1 be the areas of the sections of 

the solid then 


V=$h[(Ai + An+1) +2 (Az + As + Ar +...) 
+4(4, +4, + Ag+... .)] 
10.7. CENTRE OF GRAVITY 


The coordinates (x,y) of the centre of gravity of an area are 


p= 4 {” pa |” .: 
a ay ely OO Maced 


10.8. MOMENTS OF INERTIA 


Let I = Mk? be the moment of inertia of a body of mass M 
about a given axis; then k, which has the dimension of a length, 
is called its radius of gyration about that axis. 

Let Ioz, Joy be the moments of inertia of the plane area A of 
uniform density o about Ox, Oy, then 


Z2 vt 
Pe ) oty® dx Toy = | ox*y dx 
ite TL 


Theorem of Perpendicular Axes (applicable to plane laminae) 
Let Ox, Oy be a pair of perpendicular axes lying in the plane of 
an area A and let Oz be an axis perpendicular to both, then 


loz ici Tor + Toy 
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Theorem of Parallel Axes Let OO’, GG’ be a pair of parallel 
axes at a distance h apart and let GG’ pass through the centre 
of gravity of a body of mass M, then 


Too = Tag: + Mh? 


Let the solid of revolution and the surface of revolution have 
uniform densities p and o respectively then 


Solid Ing = p|” gmy* dx Toy =p a ay *(ty? + x*)dx 
Z1 Li 
Surface Ip, = 6 | 2ry*ds = 6 | = 2ry>/(1 + y’*)dx 
Z1 
Loy = o|™ 2ry(sy? + x?*)/ (1 + y)dx 
1 


10.9. THEOREMS OF PAPPUS AND GULDIN 


(i) Arc AB of length s does not intersect Ox. 
G is centre of gravity of the arc. 
Surface generated by the revolution of AB about 
Ox is s X 2m7. 


(ii) Area A does not intersect Ox. G is the centre 
of gravity of A. 
Volume generated by the revolution of A about Ox 
is A X 2nV. | 


of 


Section 11 
DIFFERENTIAL EQUATIONS 


11.1. FIRST ORDER EQUATIONS 


Equation Solution 


(a) Variables Separable 


@~2 = BA Wd 
Pa, = Hy) HO LID 
£02 = 40 [ronay = [éeniae + 
Examples 
(i) Decay Function: 
dx|dt = — kx x = Ce~*! 
(ii) Tension in a belt just 
slipping over a pulley: 
dT/d0 = uT T = T, e#9 
(iii) Resisted motion under 2 


gravity. Resistance propor- v = k (1 = 67") 
tional to velocity v. 


dvldt = g — kv 
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Equation Solution 


(iv) Current i in a circuit (L,R), 
t seconds after a voltage E 
is applied to the circuit. 


di X ... Eo 
L—+Ri=E = —29 (| — er FtiL 
ead 5 aad 


(6) Solution by Integrating y = exp (—JPdx) x 
Factor [fQ exp (JPdx)dx + C] 


Y +Py=Q Integrating factor exp ({Pdx) 
x 


(P, Q are functions of x (This method may be applied 
only) to (iv) above) 


(c) Homogeneous Equations 
Pe +O =0 
dx 


(P, Q homogeneous 

functions of x,y) 

Put y = vx so that 

dyldx = v + x dvi/dx 

and equation takes the Solve as (a) above 


form x = = ¢(v) 
dx 
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11.2. SECOND ORDER EQUATIONS 


Equation Solution 
(a) d®y/dx? = f(x) y =JS[fC) ax + Ax +B 
(b) d?y/dx? = f(y) . p? = (dy/dx)? 


= 2{f(y)dy + A 


= SS — F — 
ta ae d 
i Cio et 
PaaS =ixt+B 
y 
Example 
Inverse Square Law v=dx|dt t= u(- -- 
WAS ob Verte 
dxtldt? = — pix? dx {dt = £V(2u)( -— | 
given v = 0 when x = r, 
t= 0. Now solve as 11.1 (a) 


(c) Simple Harmonic Motion 
2 
TS + rx =0 (dx|dt)? = n%(a? — x?) 
(dx/dt = 0 when x = a) 


Multiply by integrating dx|dt = + n(a? — x?)}!2 
factor 2dx/dt x = acos (nt + €) 
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Equation Solution 


Examples . 
(i) Simple Pendulum: length / 
2 
7 ie 9 = 0 6 = 4 cos {[/(g/D] t + 3} 


e. . até 0 = 0 (approx.) Period T = 27v/(I/[g) 


if 0 is ata 


(ii) Compound Pendulum: h 
and k are respectively the 
distance of the centre of 
gravity from the axisand 0 = @ cos {[/ (gh/k?)]t + e} 
the radius of gyration 
about the axis. Period T = 27/(k?/gh) 
G70. ae 
aa — +2. re sn Oe 0 


= 7+ 6 = 0 (approx.) 


if A is oe 


2 
(a) +2aD + by =0 


(a, b constants) 

Put y = Ae** 

A? + 2aA+b=0 

(auxiliary eqn.) 

Ay, Az are —a + /(a? — b) Continued on p. 56 
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Equation Solution 


Case I. a® > b y =e" {A exp [V/ (a? — b)x] 
1, 42, real and unequal. + Bexp [—+/(a? — b)x]} 
Ifb =0, y= A + Be~? 
Ifa =0,b = —k? <0, 
y = Reoshkx + S sinh kx 


Case II. a? = b. 

Ay, Ao, real and equal. y=(R+ Sxe-* 
Case Ill. a? — b = —p? <0. 

Ay, 42, complex. y = Re~* sin (px + 8) 


Conjugates —a + ip. 


11.3. LINEAR DIFFERENTIAL EQUATIONS OF ORDER 
n WITH CONSTANT COEFFICIENTS 


which may be written 
(a,D" + @,-1 D1 +...+a@aD+a)y=0o0rf(D)y=0 @ 


where D = d/dx. 
Put y = Ae**, Auxiliary equation: f(A) = 0 which determines 
n values of A so that 


n 
y= p> A,e*;* 
=] 


; 
Example: Whirling of shafts. £2 — nty = 0 
xample: Whirling Of s a ee 


Auxiliary equation: A* — n* =0 
y = Ae™ + Be-™ + Reosnx + S sin nx 


56 


11.4. PARTICULAR INTEGRALS 
The general solution of 
f(D)y = $@) (ii) 


is y = u + v where u is the general solution of (i) the comple- 
gaue function and v is a particular integral of (ii) given by 


= 3B D) —— $(x) 


y = Complementary function + Particular integral 


Particular integrals may be found by the use of the following 
operational] results. 


vt ee az I 
f(D)e* = fae 7D Be “a a nd 
f(D) Soa} (ax +B) = f—a SM) (ax + 6) 
1 a ee es am eos) arse) 
f(D) ©0s f(—a’) 5 Cos 


f(D)le*V(x)] = e* f(D + a) V(x) 


11.5. FURTHER SECOND ORDER EQUATIONS 


Equation Solution 
da?x 
(a) a ~ + nx =o x = asin (nt + €) + c/n? 
(c = const.) 
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Equation Solution 


atx ; : 
(6) oa +x = asin pt xX = asin (nt + e) 
(p #n) a 
P of. ne | sin pt 
The case of resonance: x = asin (nt + ¢) — vad cosnt 
p=n 2n 
2 
(c) + 2a = + bx Complementary function 
dt? dt . | 
given by 11.2 (d). 
= c sin pt Particular Integral is 
c sin (pt — €) 


1th = py" + tate] 
where tan € = 2ap/(b — p?) 


(d) Homogeneous Equations 


2 
2D + aD + by = 4(x) 


ax? dx 
Put x = e' whence 
VY 
ax @ 
x? dy? = d*y — dy 
Ge AR ae 
d*y dy 
pr I—+b 
pS +(@-—1) an + by 
= ¢(e’) Solution as in 11.4. 
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Section 12 
STATISTICS 


12.1. DEFINITIONS 


Statistical Data A set of measured (numerical) values of x 
called the variable or variate (continuous or discrete). 


Population The totality of all individual values of the variable. 


Sample A finite set of observed values drawn from the popu- 
lation. 


Frequency Distribution Let the range of the variable be divided 
into n class intervals. Let N be the total frequency of individual 
values and let n, be the frequency in the ith class, then 


nh n 
N=2n,_ - Relative frequency f, = - afi 1 
i=l N i=} 


Statistical data may be displayed by (a) frequency and relative 
frequency curves, (6) histograms, (c) frequency polygons. 


Mode Value of x for which the frequency is a maximum. 


Median Value of x at which the population is divided into 
two numerically equal parts. 


Quartiles Three values of x which divide the population into 
four equal groups. 

Percentiles Values of x which divide the population into 100 
equal groups. 


Bs 


12.2. SAMPLE STATISTICS 


Random sample of size N. Number of class intervals n. 


n n 
rth moment about the origin = m,’ = 2. mx" = > fixie 
N j=1 i=1 
Mean X = ist moment about the origin = m,’ = & fx, 
1 
rth moment about the mean = m, = rv ym, (x; — X)" 


” Lf; (x; — ¥)° 
my, = 0 and Me oe s? =m Lf (x: — X)? = 2X fix,? — x? 
mz = Sample variance s = Sample standard deviation 


12.3. POPULATION PARAMETERS 


Similar measures for the whole population are 


44 = population mean o? = population variance 
o = standard deviation. 


For a continuous variable x having a relative frequency or a 
probability density function (p.d.f) y = f(x), 


+ 
rth moment about the origin = y,’ = - x’ f (x)dx 
+ 00 
Mean = Ist moment about the origin = 4,’ = w= xf (x)dx 
- 0 


+ 00 

rth moment about the mean = yw, = | (x — mw)" f(x)dx 
—- 00 

Population variance = fg = 0? = 2’ — pb? 


= | °° (= w? fede 


12.4. ELEMENTARY PROBABILITY THEORY 


An event comprises an outcome or a group of outcomes of an 
experiment. Outcomes may be represented by sample points in 
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appropriate coordinates. Sample space comprises the set of all 
such sample points. 

Let p; @ = 1, 2, . . ., m) be the probabilities assigned to n 
sample points, then & p, = 1. 


i 
Probability of Event A = P(A) = Xp, 


A 
the sum being taken over the probabilities of the set A of those 
sample points associated with the occurrence of the event A. 


Compound Probability Git 
P(A U B) = probability of the occurrence of the event A or B 
or both. 


P(A ( B) = probability of the simultaneous occurrence of 
events A and B. 


P(A U B) = P(A) + P(B) — P(A 1 B) 


Events A and B are mutually exclusive if the occurrence of one 
precludes that of the other, i.e. A © B = ¢. Therefore 


P(A U B) = P(A) + P(B) 


If Ai, Ao, As, . . . are mutually exclusive events, the probabi- 
lity of the occurrence of at least one is ): P(A)). 
i 


Conditional Probability 


The probability of the occurrence of event B subject to the 
condition that event A has occurred is a conditional probability, 
denoted by P(B/A). 


P(B/A) = P(A A B)/P(A) or 
P(A CB) = P(A) P(B/A) = P(B) P(A/B) 


If events A and B are independent, P(B/A) = P(B) so that 
P(A A B) = P(A) P(B) 
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If the events Ai, Az, As, . . . are independent, the probability 
of their simultaneous occurrence is 


P(Ax) P(A2) P(As). . . = I P(A) 


Bayes’ Theorem 

Suppose an event B can arise from any one of the initial or 
“causal” events Ai, Ao, . . ., An. The probability that a par- 
ticular causal event A, has resulted in the event B is given by 


P(A,|B) = P(A,) P(BIA,)| % P(A) PCBLAD 


f(,y) is the Joint Probability Density Function of x and y 
if f(x,y) is the probability of the simultaneous occurrence of 
values x and y. 


+0 [+o 
ip oO f@,y) dx dy =1 
If x and y are independent variables 
fy) = fi) fa) 


Product Moments of a Bivariate Distribution 
With respect to x = y = 0, 


i + [+00 
a: J - © J oo TEM ax dy 
With respect to the means ju, and jm, of x and y, 
+00 [+0 
Hrs = i 2 ii ve ere Mz) "(y or My) sf (x, y)dx dy 
If x and y are dependent (correlated) variables 


Covariance = cov (x,y) = Mi 4 0. 
Correlation coefficient = p(x, y) = “11/020, 


o,, 0, being the standard deviations of x and y. 
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12.5. VARIANCE OF A LINEAR COMBINATION OF 
RANDOM VARIABLES 


If z = Lia,x,, the mean Z = pw, = Da. 
i i 
var z= 0,7 = Yao? + Ldia,apyo,o; 
i {JZ 


if the variables x, are correlated. 
= La,o,? if the variables x, are independent. 
i 


12.6. DISTRIBUTION OF SAMPLE MEANS 

Random sample, size n — x,(@i= 1, 2,. . ., n) 

Sample mean = % = : 2axe 

Mean of sample means = we = 

Variance of sample means = (ox)? = 0?/n 

12.7. THEORETICAL FREQUENCY DISTRIBUTIONS 
(a) Binomial Distribution (Discrete Variable) 

Probability of x “successes” in x trials = f(x) = (3) gt 
where p = probability of “success’ and g = (1 — p). 
Mean = np. Variance = npq. 


qq hp = :. (7) pgs =» f@) 
r=0 zr=0 


In N samples, each of n trials, the frequencies of 0, 1, 2 
nm successes are terms of the expansion of N (g + p)". 


The proportion of successes in n trials = z = x/n. 


p.df. of z = a pregna-®) 


e e 
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Mean of z = (wp) = 


Variance of z = 5 (npq) = pq|n. 


(b) Poisson Distribution (Discrete Variable) 


Derived as an approximation of the Binomial Distribution. 
Let m = np be the mean number of successes in 7 trials and let 1 
become large and p = m/n small. 


Probability of x successes = e~™ m'*/x! 
Mean = Variance = m 


(c) Normal Distribution (Continuous Variable) 
Denoted by N[p, ss Mean uw. St. dev. o. 

x= 2 
bat ie = m7 a5? (pea 
(—-7o<x< +o) 
Standard Normal Variable = u = (x — p)/o 
Normal Probability Integral for the calculation of cumulative 
normal probabilities 


exp (—4u?) du is tabulated. 
Teas Jo 8? (40 
12.8. ELEMENTARY SAMPLING 


(a) If samples are drawn from a normal population, the sample 
means are normally distributed. 

(b) Let x1 be N[f1, o1] and x2 be N[f2, o2]—a pair of inde- 
pendent variables—then a,x; + d2X2 is 


Nla@ipfir + d2fle, V/(a12017 + ag?027)] 
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(c) Central Limit Theorem Sample size n drawn from a popu- 
lation of mean y, st. dev. o. 


Sample mean X and y = ( — u)(V/n)/o 
p.d.f. of y — N[0, 1] as n— 0 


(d) Confidence Limits for the mean y of distribution N [u, o] 
from a sample, size n, mean X. 


X¥+1960//n (95%) ¥A2.58 a//n (99%) 


12.9. SMALL SAMPLE DISTRIBUTIONS 


(a) Small sample yi = 1, 2, .. . ., 2) drawn from population 
N[O, 1]. 


n 
u= 4? = & y? has a x? distribution f(y?) with n degrees of 
i=l 
freedom (d.0.f.) given by 


Bt ms (y?)i"- 1e- 4x2 
I) = f= ~ Tanz 


The cumulative x? distribution i f) du has been tabulated. 


(b) Confidence Limits for the variance o? of a normal variable. 
Sample, size n, variance s?. (100 — «)% confidence limits for 
o? are 


ns?/4,7 > 02 > ns?/4.? 
where P(x? < 417) and P(y? > yo?) are each $a °%. 


(c) F-distribution If u, v are independently distributed y? vari- 
ables with 71, v2 d.o.f. then if F = u/v, + v/ve, the p.d.f. p(F) is 


F 
known and the cumulative distribution function P(F)dF has 


been tabulated for various d.o.f. 
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(d) F-ratio Test 


To test, on the basis of variance, whether two samples, sizes n,, 
Nz, Variances $17, 52”, are likely to have been drawn from the 
same population of variance o*. 


U = N41812/07, v = NeS22/o7 have x? distributions with 
(n, — 1), (n2 — 1) dof. 


M1837 / NeS2” A 
Fa —— | ———_ = 6,"/6," 
(ny — 1)/ (2 — 1) oo 


where 6,2, 62? are unbiased estimates of o?, on the null hypothesis 
that the samples are drawn from the same population. 


(e) f-distribution 
If u is N[O, 1] and v has a x? distribution with » d.o.f., the 
probability distribution of t = uv/»//v has been tabulated. 
Sample size n drawn from N[p, ol]. 
Taking u=(% —)Vn/o [see 12.8 (©)] and v = ns?/o? 
having a 7? distribution with n — 1 do. 
t= (% — w)V(n — 1)/s has a t distribution with n — 1 d.o.f. 


(f) Confidence Limits of the Mean 


if sie. acacia Ry CON ee la j2 [1001 —a %] 


s 
V(n — 1) V(n — 1) 
where P(t > ta)2) is «/2. 


(g) t-test: Difference of Two Sample Means: Is there any essen- 
tial difference between two similar samples of data? Suppose 
x,y are distributed N[:, o;] and N[u,, oy) respectively. Two 
random samples, sizes n,,/,, means x,y. 

On the assumption that the samples are drawn from the same 
population, x and y are distributed with common variance 
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o = 0; = Gy [justify the assumption by use of F-test, 12.9(d)] 
and fd, = py. 


(a Se ( j wtiea ) with (n, + m, — 2) dof. 
Nz + Ny 


6 
and 6 is the unbiased estimate of o from pooled data. 


12.10. GOODNESS OF FIT 


Population divided into k mutually exclusive classes. f;, m, are 
respectively the observed and expected frequencies in class i. 


3 Cy ih)” 


fae ores 


has a x? distribution with (k — 1) d.of. 


12.11 LINEAR REGRESSION 


Bivariate distribution. Sample size N. Sample points: (2, y,). 
Frequency /f, and N = 2 f; 


(a) Least Squares Regression Line y = a + bx 
a ‘ Nts ae SS 
a =f — bx where x = Hy fed = yh 


NX fixiy: — DL fixiaf; 
and 6 = ———+____—_—+—— 
NSfx® — fan? 


(b) Unbiased Estimate of Population Correlation Ratio p 


is given by 
yo gt Is phage x fix: — X) (yi — J) 
Mn eM CR ae 
where 5,’, Sy’, Szy’ are sample st.devs. and covariance respectively. 


(c) Significance of r On the hypothesis that p = 0, 
r/[(n — 2)/( — r?)] has a t-distribution with (n — 2) d.o.f. 
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If the sample value of f is significant, the hypothesis is rejected 
and a degree of correlation is accepted. 


(dq) Rank Correlation (Spearman) Items ranked and assigned 
ordinal numbers in respect of some characteristic (usually 
magnitude). 

Values of x have ranks u, and values of y have ranks v. 


62d,? 


Spearman’s coefficient r= 1 — 
P (n®? — n) 


where d, = u,; — 1. 
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Section 13 
NUMERICAL ANALYSIS 


13.1. ROUNDING OFF ERRORS 


Let A, be the rounded-off measured value of a quantity whose 
exact value is A,. 
(Maximum rounding error = 3 unit of the last decimal 
place in A,) 
Absolute error = a = A, — A,. 
Absolute error modulus of A, = {a}. 
Relative error = a/A, 


Let a1, a2, ag be the absolute errors of three measured values 
(A,)1, (Ax)a, (Ap)s respectively, then if 


(A,)3 = (A;y)1 & (Aro las] < Jai] + Jae] 


(A re = (A,)1 + (A;)o jae hes. ane ae 
or (A,)1 X (A,)2 ie rs (A;)3 (Ar)o 


(Ayo = {((Ar 1}? saa |p M 


cm 


13.2. ERRORS IN FUNCTIONS 


(a) If u = f(x) and x is subject to error, then 


éu = Dx 
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where du, 6x are the absolute errors in u, x respectively. 
(b) If u = f(x, y,z.. .) then 


of of of 
du=—_6d eG ee ae 
a on 


13.3. TABULAR VALUES OF f(x); DIFFERENCE 


NOTATION 
Constant tabular interval = A 
Value of x Value of f(x) Notation 
X-~2 feo — 2h) f-2 
X14 f (x iar h) sf -1 
Xo yt (x o) f 0 
x1 sf ro +h) f 1 


X2 fo + Zit) fa 


Forward Difference Operator A 

First difference: Afn = Af (xn) = f(xn + 4) — fxn) 
Second difference: A(Af,) = A?f, = Afna1 — Afn 
rth difference: Af = AP *faar — Atm, 


13.4. LINEAR INTERPOLATION 
If x, = Xo + ph, then f, = fo + p Afo 
13.5. SHIFT OPERATOR E 


Ef f(x)} =f + 4) EX f@} = ELE f@)} 
= E{ f(x + h)} = f& + 2h) 
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In general 

E*{ f(x)} = f(x + ph) (p +ve, —ve or fractional) 
A more concise notation for this is 

E’fn = fn+p 
In particular, Ef_1 = fo . EO aren 

Ee fo = fra = f(x + 44) 


13.6. RELATIONSHIP BETWEEN A AND E 
Afn = foi — fn ng (1 + A)fn = fn+1 = Efn 
* E=1+A and E£'=(1 +A) 
(r +ve integer) 
13.7. GREGORY-NEWTON FORWARD DIFFERENCE 
INTERPOLATION FORMULA 
Value of f(x) required at x = x9 + ph. 
fr = E*fo = (1 + Ao 


—1 —1)\(p—2 
nam | f + pAfo + Pr ) arp, PP ww ) Ae, 


Nee 
13.8. BACKWARD DIFFERENCE OPERATOR V 


Defined by Vf(Qn) = fn) — fn — A) or Vin = fa — fra-1 
and Vf, = V'~*f, — V'~*fn-1 

E-*f;, = fn=1 = fn ay Vin =(1-— V)fn 

E-*=(1—V) or E=(1-—V)7? 


ry 


13.9. GREGORY-NEWTON BACKWARD DIFFERENCE 
- INTERPOLATION FORMULA 


fo = Erfy = (1 — V)-*fo = {fo + pV fo Pere )) yap, 
Pip + DY +2) a 
+e+ 00+ Dame, | 


13.10. CENTRAL DIFFERENCE OPERATOR 8 


6 is defined by Ofn+i2 =fn+1 — fn for all integral values of 
n. The operator 6 may be applied to any column of the difference 
table so that 


Of = Ofn+1/2 hiss Ofn=112 

O'fn+1/2 ol O7fn+1 ean O*fn etc. 
In particular, 6?"*? f,,2 = 6*f, — 6"fp (r= 1,2,3,...) i) 
Ofn+12 = fnsi — fn = EB fasije2 — Eo? fasiye 

= (E12 — Brot) foi 

Therefore 6 = E1/?2 — E-3/2 

E¥26=(E—1)=A 62 =A2/E= A211 + A) 
13.11. CENTRAL DIFFERENCE INTERPOLATION 

FORMULAE 


These formulae use information placed symmetrically with 
respect to the interpolation interval. 


(i) Bessel f, = 


fo + P fara + PO (%fo + d7f1) +L =P =D 599, 


$ Pt OP Ne (aif + Of) +. ni 
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(ii) Everett 


Using (i) of 13.10, the odd order differences may be eliminated 
from the Bessel formula above giving 


fo = (1 — pif + pfs + 2L— PEP) 509, 
Pm il horn ow + 1) fF, +e 
[Further coefficients will be found on p. 56 of Interpolation and 
Allied Tables (I.A.T.), H.M.S.O.] 
(iii) Stirling 


fo = fo + Si(Of-1)2 + Oftj2) + S267fo + Ss (O?f—1)2 + O*fi)2) 
+ S.0*fo t+... 


The coefficients S;, So,. . . will be found on p. 58 of £.A.T. 


13.12. LAGRANGIAN INTERPOLATION 


f(x) takes" the values «fo, fs. sss, ig TOT Wee XG HA Ps en 
respectively. It is required to find a polynomial which takes 
these (n + 1) functional values for the same values of x. The 
required interpolate is found by calculating the value of this 
polynomial for the given value of x. 


Lagrangian Interpolation Coefficients 


Lx) = 
(% — Ror Ri) i ee em a). A me Xa 
(x4 — XoNXy — X1) oo gm Kya Km p41). 2. KE Xn) 


defines the (n + 1) coefficients for i = 0, 1,2,...n. 
L(x) are polynomials of degree n. 
L(x;) = 04 ff) and L,(x,) = 1. 
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The Lagrangian Polynomial is defined by 
nm 
| L(x) = UL) f,_ which is of degree < n 
i=0 


Since L(x,) = f, for i=0, 1, 2,. . ., m, L(x) is the required 
polynomial. If f(x) is of degree <n, it may be shown that 
f(@) = LQ). If f@) =c, 


LX) = ¢= > L(x) .¢ Ss L(x) =1 (check) 
i=0 i= 0 


13.13. AITKEN’S METHOD FOR INTERPOLATION 
WITH UNEQUAL INTERVALS 


f@) has values fo, fi, fj at x = x0, X1, xy respectively (x1 — Xo 
# X; — X1, in general). Determine the interpolate at x. 


(x id Xo) fo 
(x ¥3) fy 
(xs — Xo) 

(x — Xi) foil) 

(x — x). . ie 


Second tion = | 
econd approximation = fo1; (x) EEE, 


First approximation = fo,(x) = 


13.14. NUMERICAL INTEGRATION 
If x = xo + ph hs f(x)dx =h jhe dp 
nae) 


(i) Forward Difference Integration Formula (derived from 
Gregory—Newton formula) 


[Presa = [nfo +2 a f+ 20 


a asf 


oe 
+ apt. ..| 
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and, in particular, 

fr ‘foods =h| fot Wfo—LAaye+s A%+...| 
Simpson's Rule (an even number n equal subdivisions of width h) 
[2 food = AM fe + fa + fa that. t faa) 

TA. + fs Se fai) 

Estimate of error = — = 9 tA‘fo +Atfiit+...+ A*f,-2} 
Trapezoidal Rule (n equal subdivisions of width h) 

[2 foode = ait fo + fit fat. + fad +f 

Estimate of error = 

“ = (uefa — 40 fo) ++ a= 55 (fa HO°fo) 


where 4 is the averaging operator defined by 


Mfp _ 4 (fp+1/2 + fp-112); 
LO" f, = 3(0"fp +12 + O"fp-1:2) (.A.T. p. 54) 


(ii) Central Difference Integration Formulae 
(a) Derived from Bessel’s interpolation formula. 


[P feddx =n [Ae +4) - ae + 8f,) 
+o ; = (Of) + Of) +.. | 
=h{nfia—s 5 MP faia + Oise + ‘ : 
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(6) Derived from Stirling’s i Ne formula. 
z4 
= = 0*fp — —~ 64 ae " 
[P feoa 2h (fo +2 *fo “a o+ 
13.15. NUMERICAL DIFFERENTIATION 
/ / d, l / 
Hx=xotph f@M=f(pZ=-f—) 


Se —— 


(i) Central Difference Formulae 
(a) Derived from Bessel’s Interpolation Formula. 


bie ¢: =< {fiat 4(2p om 1)(0°fo ~ 07/1) + 6(3p? AN 3p + 3) 0% fj j2 
+55 Op — Dip =p — NOY + Of) +. . 
At half-way points f’(xo + 3h) 
1 1 3 
ah {ofan — 7A Of 12 a Fag ° fue aa oe i 
(5) Derived from Stirling’s eae Formula. 
f'0) = 2 {HOf- aia + Sus) + po%fo + 4 Gp? — 16-1 
+ Pfura) + 5 pOp? — DH fot.. 4 
At a tabular point xo, p = 0 so that | 
f (Xo) =; {45f- 11 i Of j2) rT 5 (8°f-119 a O° fi ;2) 
_ sh +O Fie Tr <5 : 


; {udfo — 2 oe HO'fo + = 5 Hfo + an 
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(ii) Forward and Backward Difference Formulae (derived from 
the Gregory~Newton Formulae) 


(a) Forward Difference formula. 


(0) = 7 |Afo + 4Op — WA + AGp? — Gp + DAY, 


+ oe 


For the tabular point xo, put p = 0 so that the Forward Differ- 
ence formula becomes 


Fo) =1{Afy— Ay + MY —1aY, +2 A%—.. | 


(b) For a point near the end of the table, use the Backward 
Difference formula 


fia) = (Vf + AVS + IVS + AV +205 +... 
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